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The charged particle stopping power in a highly ionized and weakly to moderately 
coupled plasma has been calculated exactly to leading and next-to-leading accuracy 
in the plasma density by Brown, Preston, and Singleton (BPS). Since the calcula- 
tional techniques of BPS might be unfamiliar to some, and since the same methodol- 
ogy can also be used for other energy transport phenomena, we will review the main 
ideas behind the calculation. BPS used their stopping power calculation to derive a 
Fokker-Planck equation, also accurate to leading and next-to-leading orders, and we 
will also review this. We use this Fokker-Planck equation to compute the electron- 
ion energy partitioning of a charged particle traversing a plasma. The motivation 
for this application is ignition for inertial confinement fusion — more energy deliv- 
ered to the ions means a better chance of ignition, and conversely. It is therefore 
important to calculate the fractional energy loss to electrons and ions as accurately 
as possible, as this could have implications for the Laser Megajoule (LMJ) facility in 
France and the National Ignition Facility (NIF) in the United States. The traditional 
method by which one calculates the electron-ion energy splitting of a charged par- 
ticle traversing a plasma involves integrating the stopping power dE/dx. However, 
as the charged particle slows down and becomes thermalized into the background 
plasma, this method of calculating the electron-ion energy splitting breaks down. 
As a result, the method suffers a systematic error of order T/Eq, where T is the 
plasma temperature and Eq is the initial energy of the charged particle. In the case 
of DT fusion, for example, this can lead to uncertainties as high as 10% or so. The 
formalism presented here is designed to account for the thermalization process, and 
in contrast, it provides results that are near-exact. 



I. INTRODUCTION 

When a fast charged particle with initial energy Eq traverses a plasma, it looses energy at 
a rate dE/dx per unit of distance, and it comes into thermal equilibrium after depositing its 
energy into the electrons and ions that make up the plasma. Using the formalism of Brown, 
Preston, and Singleton (BPS) |l|, which provides a means of regulating the kinetic equations 
at short and long distances in a consistent manner while treating quantum mechanical effects 
exactly, we shall calculate this energy splitting in a uniform plasma to leading and next-to- 
leading order in the plasma coupling. Our motivation is the thermonuclear burn of deuterium 
and tritium in inertial confined fusion experiments where a fast a particle is born with an 



energy Eq = 3.54 MeV. The more efficiently the ions are heated, the easier it will be to 
initiate the bootstrap heating process that triggers ignition and burn. 

In this paper we shall always work with a plasma whose components have a single tem- 
perature T (measured in energy units). We are currently generalizingour formalism to a 
plasma in which the electrons and ions are at different temperatures [2[] . The energy par- 
titioning is usually computed within the context of a fast charged particle traversing the 
plasma until coming to a complete stop, in which case the energy partition into ions and 
electrons is given by 

E <=L dE iik (L1) 

and 

E [ E °dE dEe / dx (1 2) 

e J dE/dx 

Here dE l /dx and dE e /dx are the stopping power contributions from the ions and electrons, 

dE _ dE I + dE e 
dx dx dx ' 

and thus E t + E e = Eq. This is only an approximate description because the fast charged 
particle does not simply come to rest within the plasma, but rather it becomes thermalized 
at the temperature T. One should not extend the integrals in Eqs. (11 .ip and (11.21) down to 
zero, but rather a lower limit E m \ n on the order of the thermal plasma energy, -Emm ~ T. 
Consequently the systematic error in the above calculation of E l and E e is of order T/E , 
and as we shall see, the correct electron-ion energy partition relation reads 

E l + E e + ^T = E . (1.4) 

Before examining the energy splitting, we need to discuss the calculational framework 
within which it appears. As we shall see, the correct expression for the energy splitting 
arises from a Fokker-Planck equation derived in BPS; however, before jumping into details, 
it will be useful to briefly review some salient features of the stopping power calculation. 



II. THE BPS FORMALISM 



Calculating Coulomb energy exchange processes in a hot plasma is notoriously difficult 
because of the subtleties of the Coulomb interaction, which produce logarithmic divergences 
at both long and short distance scales. This problem was first spelled out and solved 
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to leading order by Landau and then Spitzer in the context of electron-ion temperature 



[flfl. 



equilibration, and later by Corman et al. for the charged particle stopping power 
Since the divergences are only logarithmic, one introduces ad hoc short and long distance 
cutoffs b min and b max , and the rate of energy loss of some process (such as temperature 
equilibration or stopping power) can be cast in the form 

Coulomb Logarithm 

The prefactor K is easy to compute exactly. The logarithmic term, conventionally called the 
Coulomb logarithm, can only be approximated within the above scheme. The long distance 
scale 6 max is set by the relevant Debye screening length, while the short distance scale 6 m i n 
is determined by either the Landau length or the thermal de Broglie wave length (or some 
interpolation between them). As such, this method suffers a systematic uncertainty in the 
argument of the Coulomb logarithm. 1 In the language of perturbation theory, Eq. (12. ip 
is accurate to leading order in the plasma coupling constant. This accuracy was extended 
to subleading order by BPS [l| which performed a first principles controlled calculation, 
including the exact terms under the logarithm and a rigorous treatment of the quantum to 
classical transition. 

Reference jg] discusses at length the manner by which one can expand thermodynamic 
quantities such as dS/dt as a perturbation series in powers of a small parameter g, the 
plasma coupling defined by 

9 = ^, (2-2) 

where k is the Debye wave number. This is just the ratio of the potential energy of two 
electrons a Debye distance apart to the thermal kinetic energy of the plasma, and it is 
related to the usual plasma parameter by g oc T 3//2 . Quantities expand in integer powers of 
g, except for possible In g terms, and the rate of energy exchange for the stopping power 
takes the form 

^ = - Ag^ng + B£ + 0(g 3 ) = A g 2 In A coul + 0(g 3 ) . (2.3) 

LO NLO 



1 The constant under the logarithm sometimes varies by an order of magnitude from paper to paper within 
the literature, depending upon the choices of 6 m i n and 6 max - 



3 



We have indicated the leading order (LO) and the next-to-leading order (NLO) terms in 
the ^-expansion. Here lnA coul = — ln{Cg}, with C defined by B = —A InC, and therefore 
knowing the next-to-leading order coefficient B is equivalent to knowing the exact coefficient 
C under the logarithm. To get a feel for the numbers, at the center of the sun g = 0.04, and 
the error term in Eq. (12. 3p is consequently small. 

The problem with directly calculating A and B is that the kinetic equations diverge and 
must be regularized in the appropriate manner. Furthermore, to find the coefficient under 
the logarithm, this regularization procedure must preserve the delicate balance between 
the long and short distance physics. Indeed, the BPS calculation 1] includes both short 
distance physics and dynamic collective long distance physics, joined together exactly and 
unambiguously (and this is the rub), systematized by a power series expansion in the plasma 
coupling constant g. The coefficients are also calculated to all orders in the dimensionless 
quantum two-body scattering parameter rj, thereby providing an exact interpolation between 
the extrme classical and quantum regimes. 

The rigorous starting point is the BBGKY hierarchy (or its quantum generalization), 
which is finite and well defined and does not suffer from the aforementioned divergences. 
One must of course truncate this vast number of equations to something manageable, such 
as the Boltzmann or Lenard-Balescu equations, and it is this truncation process that renders 
the various three dimensional integrals divergent. However, as shown in Ref. (7|, in v spatial 
dimensions these divergences become simple poles of the form 1/(^—3). In spatial dimensions 
v > 3 the BBGKY hierarchy reduces to the Boltzmann equation (BE) to leading order in g 
(the BE is finite and does not have the usual long distance divergence for v > 3). Calculating 
the rate of energy loss using the i/-dimensional BE gives a result of the form 



The "greater-than" superscript is to remind us that the calculation has been performed in 
dimensions v > 3. In a similar manner, to leading order in g the BBGKY hierarchy reduces 
to the Lenard-Balescu equations (LBE) for v < 3 (the LBE is finite and does not suffer from 
short distance divergences when v < 3). A calculation of the energy rate with the LBE gives 
a form 




LO in g when v > 3 . 



(2.4) 




LO in g when v < 3 . 



(2.5) 
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Note that both rates are of order g 2 in three dimensions, and they both suffer from a divergent 
simple pole. The coefficients H(v) and G(v) can be expanded in powers of e = v — 3, with 

H(u) = -A + eH 1 + 0{e 2 ) and G{v) = -A + eG 1 + 0{e 2 ) . (2.6) 

The leading terms must be equal, H{3) = G(3) = —A This arises from the calculation 
itself and is not imposed by hand, and it makes the short- and long-distance poles cancel, 
thereby giving a finite result. 

Since the rates d£ > /dt and d£ < /dt were calculated in mutually exclusive dimensional 
regimes, one might think that they cannot be compared. However (and this is perhaps the 
most crucial step in the method, and certainly the most subtle) , we can analytically continue 
the quantity dS < /dx to dimensional values v > 3, after which we can directly compare the 
rates (12.41) and (12. 5p in a common dimension u, and the limit v — > 3 may then be taken. 
Upon writing the ^-dependence of Eq. (12.51) as g 2+ ( u ~ 3 \ when v > 3 we see that the rate 
(12. 5p is indeed higher order in g than Eq. ( 12. 4ft since e = v — 3 > 0: 

—— = -G{u) + 0(u - 3) : NLO in g when v > 3 . (2.7) 

The individual pole-terms in Eqs. (12.41) and (12. 7p will cancel giving a finite result when the 
leading and next-to-leading order terms are added. Summing terms ( 12. 4ft and (12.71) . using 
the relation g e = expjeln g} = 1 + elng + 0(e 2 ), and taking the e — > limit gives 

^ = -Ag 2 \ng + Bg 2 + 0(g*), (2.8) 

with B = H\ — G\. This is in agreement with Eq. ( 12. 3p . In this way, BPS has calculated the 
charged particle stopping power accurate to leading order and next-to-leading order in g. 

BPS also derived a Fokker-Planck equation accurate to leading and next-to-leading order 
in the plasma coupling [1]. Denoting the phase space density for the dilute collection of 
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/(r,P,t), (2.9) 



charged particles by /(r, p, £), this equation reads 



where the sum runs over the plasma components b, (3 = 1/T, the vector v = p/m is the 
velocity of a particle with momentum p, and the summation convention is used for repeated 
vector indices. The symmetric tensor C^ 1 has longitudinal and transverse components, 

Cf (p) = ME) ^ + B h {E) \ (6 H - v k v e ) , (2.10) 



where v = |v| and v = v/v, and E = |mt> 2 . We denote the sum of the ion components 
of the ^.-coefficients by A Y = Yli^i and the electron component by A e , with A = A Y + A e . 
Expressions for the Ab can be found in BPS [l]. With our conventions, the number and 



kinetic energy densities of the charged particles is given by 

d 3 p 



n(r, t) 



f(r,P,t) 



and 



S(r,t) 



(2.11) 



(2.12) 



(2tt^) 3 2m 

We can derive a relation between the stopping power and the ^.-coefficients in the fol- 
lowing manner. For a single particle at r p moving with velocity v p , the distribution function 
takes the form f(r, p, t) = (2nh) 3 5(r — r p )5(p — p p ), and the Fokker-Planck equation gives 
the particle's rate of energy loss as 



dE 

dx 



dE 
~dt 



a i 9 



<ct\p P 



(2.13) 



Upon substituting the decomposition (12.101) for the scattering tensor and dropping the pro- 
jectile subscript p, the contribution from species b appears as 



dE b 
dx 



d 



(3mv dv e 



Ab(E) . 



(2.14) 



As v gets large (f3mv 2 ^> 1) note that dE^/dx — > Ab 



III. FORMULATION OF THE PROBLEM 



With this background in hand, we turn now to the energy splitting problem. Rather 
than tracking an individual charged particle slowing down in the plasma, it is much simpler 
to examine a homogeneous and isotropic source of charged particles of a single energy E . 
The distribution function will therefore depend only upon the energy and time, / = f(E, t). 
Furthermore, the homogeneity and isotropy conditions will greatly simplify the form of the 
Fokker-Planck equation: only the transverse coefficients Ab will enter the diffusion kernel 
on the right-hand side of Eq. ( 12. 9ft . while the convective term on the left-hand side will of 
course vanish. Thus, we now consider the inhomogeneous Fokker-Planck equation with a 
source, which, after some algebra reads 



§f-—-ipME)E 
at mv oE 



d 

1 + T — 

8E 



f(E,t) = 6(E-E ) s(t) . 



(3.1) 



We find [2i] that the relevant solution is of the form 

f(E,t)=n(t)Afe-? E + f(E) 



(3.2) 



We proceed now to motivate the structure of the solution (13.21) and the nature of the time in- 
dependent function f(E). The particles will eventually thermalize to a Maxwell-Boltzmann 
distribution, and the first term of Eq. (13.21) merely represents the buildup and subsequent 
thermalization of the particles produced by the source. The normalization factor M is chosen 
so that n(t) is the number density of the produced particles once they have thermalized into 
the Maxwell-Boltzmann distribution oc e~ /3E . The time-independent piece f{E) describes 
the steady state of nonthermal particles losing energy to the plasma, i. e. particles cascading 
down "energy bins" from the initial energy Eq to the final thermal energy. The situation 
described here can be pictured as the flow of water over a rocky waterfall that slows the 
motion of the water as it descends. The initial rate of flow of the river corresponds to the 
rate of produces particles; the height of the waterfall to the initial energy Eq. The energy 
dissipated in the fall corresponds to the energy lost to the ions and electrons and is deter- 
mined by f(E). The final flow into a horizontal lake corresponds to the build up of the 
particles into their final thermal equilibrium state. This is illustrated in Fig. [TJ 

We now specialize to the relevant case in which the source s(t) slowly turns on and attains 
a constant value. In this case, particles become produced at a constant rate per unit volume 




n(t) 



FIG. 1: The waterfall analogy. The small green rocks represent the plasma electrons while the 
larger red rocks are the plasma ions, with the flowing 'water' in blue representing the evolution of 
the produced charged particles (a particles for DT fusion). As 'water' falls down the electron- ion 
slope at a constant rate determined by f(E), energy is deposited into electrons and ions. At the 
bottom of the fall is a lake into which the excess 'water' drains, representing the final thermalized 
particles, with height n(t) and Maxwell-Bolztmann distribution oc e~@ E . 
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ho, and from Eq. (12.121) the rate of change in energy density becomes 



E (t) = E h - 



{2irhf 



[Aj -\- A e 



v + 



J_d_ 

(5m dv 



RE) 



The fraction of energy lost to the ions and electrons is now identified as 



and 



Eq 



El 

E 



h E 



d 3 



P 



h E 



(2vrfr) 3 
d 3 p 



A, 



A f 



i a 

(3m dv 



1 d 



(3m dv 



RE) 



RE) . 



(3.3) 



(3.4) 



(3.5) 



(2tt^) 3 e 

In the steady state, the energy density build up of final particles is their thermal energy per 
particle times the increase in number density, E = |Tno- And so, in view of Eqs. (13.41) and 
(13. 5p . the energy balance expression (13.31) now appears as 

3 



Th = [E -Ej- E e ] h , 



(3.6) 



which gives expression (\1A\i : the original energy of a produced particle is lost to the ions 
and electrons with the remainder being the thermal energy of a free particle. 

We can simplify Eqs. (13. 4p and (13. 5p by calculating the action of the operator in square 
brackets on f(E), and one finds {2} 

" Eo dE A X (E) 



and 



El 

Eq 



Ee 

E n 



E A(E) 



Eo dE A e (E) 
~E~o A{E) 



erf(V/3#) 



A(3E 



-f3E 



7T 



(3.7) 



erf(V^E) - 



A(3E 



7T 



(3.8) 



where erf (2) is the error function. Since dE^/dx — > Ab for large energies, and the error 
function approaches unity, at high energies E we see that Eqs. (13. 7p and (13.81) approach 
the same form as the more intuitive but less accurate results (11.11) and (11.21) . The primary 
differences occur for E ~ T. 

Let us compare these near-exact results with the less precise but well known result of 
Fraley et al. (FR) [3]. Starting with a phenomenological model of the stopping power, these 
authors show that the simple rule 

1 



El 

Eq 



1 + 32 keV /T e 



(3.9) 



provided a good fit to their results. Figure [2] shows the fraction energy loss to ions for BPS 
and FR for a DT plasma with electron number density n e = 1 x 10 26 cm -3 . We see that 
FR somewhat underestimates the energy deposited to ions for temperatures up to around 



= 1 xlO 26 cm- 




50 100 150 200 

T (keV) 



FIG. 2: The fractional energy loss into ions as a function of the plasma temperature for an a particle 
in an equimolar DT plasma with initial energy Eq = 3.54 MeV. The electrons and ions have a common 
temperature T and the electron number density of the plasma is n e = 1 x 10 26 cm~ 3 . The solid line is the 
analytic calculation of this work, and the dashed line is the fit provided by Fraley et al. 
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FIG. 3: The percent change of the ion fractional energy loss between Fraley et al. and BPS for a 3.54 MeV 
a particle in an equimilar DT plasma for three densities n e = 10 24 , 10 25 , 10 26 cm -3 . 
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100 keV, and FR slightly over estimates E l for larger energies. In Fig. [3] we compare the 
percent difference between BPS and FR over a wide range of densities. 
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